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The quantum phase transition in iron-based superconductors with 'half-Dirac' node at the electron 
Fermi surface is a T = structural phase transition described in terms of nematic order. An effective 
low energy theory that describes half-Dirac nodal Fermions and their coupling to Ising nematic 
order that describes the phase transition is analyzed using renormalization group(RG) study of the 
large- Nf version of the theory. The inherent absence of Lorentz invariance of the theory leads to 
RG flow structure where the velocities vf and va at the paired half-Dirac nodes (11 and 22) in 
general flow differently under RG, implying that the nodal electron gap is deformed and the C4 
symmetry is broken, explaining the structural (orthogonal to orthorhombic) phase transition at the 
quantum critical point(QCP). The theory has Gaussian fixed point A* = 0, (va/ (vfIvf')')* = with 
stable flow lines toward it. This suggests second order quantum phase transition between Ising 
nematic ordered and disordered phases, in the presence of half-Dirac fermions in the background 
superconducting state. Interpreting the fermion-Ising nematic boson interaction as a decay process 
of nematic Ising order parameter scalar field into half-Dirac nodal fermions, we find that the theory 
behaves as systems with dynamical critical exponent z = 1. The nematic critical fluctuations lead to 
broadening of quasiparticle spectral function except at the region near the half-Dirac nodes where 
'Fermi arc-like' sharp peak of spectral weight appears. The nematic critical fluctuations do not 
modify much the anisotropy of such peak. We however found critical point A c , directly related to 
that of 4> field transition, around which the Fermi arc of spectral function collapses. 



I. INTRODUCTION 



Quantum phase transition in strongly correlated sys- 
tems such as high T c superconductors is one of the most 
active topics in condensed matter physics. In cuprates 
and the recently discovered iron-based family of high T c 
superconductors, there exists a quantum critical point 
at T = deep inside the superconducting dome that 
represents such quantum phase transition(Fig. la),b)). 
This QCP also separates tetragonal and orthorhombic 
crystal structures and thus represents structural phase 
transition at zero temperature. It has been argued 
that the orthorhombic state is described by the so-called 
Ising nematic order [1,7,8] and such structural transition 
in cuprates[2,3] and iron-based superconductors[4,5,6] 
(where d-wave symmetry was assumed) is nematic tran- 
sition. The general phase diagram of iron-based super- 
conductors illustrated in Fig. 16) shows that there is 
a tetragonal to orthorhombic structural phase transition 
at some finite temperature in the undoped case down 
to T = at a critical doping x c deep inside the dome 
where the Ising nematic order coexists with the s±-wave 
superconducting state. At T = this critical doping 
is a quantum critical point between Ising ordered state 
and Ising disordered state. The theory of quantum phase 
transition at this quantum critical point is the focus of 
this work. 




FIG.l a)Quantum phase transition in cuprates with the 
quantum critical point and quantum critical region nearby, 
b) Typical phase diagram of iron-based superconductors with 
similar quantum critical point. In both a) and b), the param- 
eter x is defined as and x c = —r and x — — r(A) as given in 
Eqs. (2) and (31). Physically, x may represent doping level, 
pressure, or other appropriate experimental quantities. 

The quantum phase transition in cuprates that re- 
lates structural phase transition with nematic order was 



first studied using renormalization group approach[ll, 
12] which showed using perturbative RG calculation at 
fixed Nt with e expansion around 3 + 1 dimensions that 
the velocity anisotropy in the nodal fermion action and 
the anisotropic coupling between nodal fermion and Ising 
nematic order leads to a fluctuation-induced first order 
phase transition, as indicated by the runaway RG flows. 
A large- Nf study of the same system but in 2 + 1 dimen- 
sions [2] however found a second order quantum phase 
transition and has finite renormalized velocity anisotropy 
as compared to Dirac-like theory such as QED 3 which 
found velocity anisotropy to be irrelevant. Another RG 
study on the same system in 2 + 1 dimensions [3] found 
vanishing velocity anisotropy ratio (^)* = as the fixed 
point. 



The coupling between nodal quasiparticles to the ne- 
matic order was argued to be the most effective driv- 
ing force of the structural transition. The presence of 
nodes and the resulting nodal quasiparticles in d-wave 
cuprates is therefore of crucial importance here. On the 
other hand, from the aspect of gap symmetry, iron-based 
family was originally thought to have isotropic s± wave 
symmetry, thus ruling out the presence of nodes. How- 
ever it was found later that the electron Fermi pocket in 
iron-based superconductors admits anisotropic gap and 
thus possible existence of nodes. 



In a related development, it has been shown re- 
cently that iron-based superconductors can have the so- 
called accidental ('zero') node( Fig. 2)) at the electron 
pocket[9,10] due to the gap anisotropy where the gap 
just touches the Fermi surface, that is, it is right at the 
onset of being gapless. In the simplest model, we can 
represent it by A(0) = A (l - cos46) (Fig. 2). Such 
accidental zero has anisotropic dispersion which in lin- 
ear in p x direction and quadratic in p y direction or vice 
versa. One can therefore interpret such zero as "half- 
Dirac" node, because it has Dirac spectrum in one di- 
rection but has parabolic dispersion in the perpendicu- 
lar direction as that for free particle. This accidental 
node is robust and persists to zero temperature where 
we can have a quantum phase transition between fully 
gapped (nodeless), zero and nodal states by tuning ap- 
propriate parameters, e.g. the coupling constants in the 
Hamiltonian. Such gapless point is accidental because 
it does not arise from or protected by symmetry. This 
kind of nodes however has recently been shown to ex- 
ist within finite regime of microscopic parameter space 
(the strength A^ of the inter-electron pockets hybridiza- 
tion) [26] and so this makes it interesting to study the 
nature of quantum phase transition in iron-based super- 
conductor compounds with such peculiar gapless point. 
This is the purpose of this work. 




FIG. 2 a)Iron-based superconductors' (reduced) Brillouin 
zone and the electron Fermi surfaces ('pockets') at (±7r,±7r) 
and hole pockets at T = (0, 0) . The electron pocket has 
anisotropic gap and its critical half-Dirac nodes paired as 11 
and 22.b)The profile of a half-Dirac node. Energy E ~ k x in 
x direction but E ~ ky in y direction. 

We summarize our main results as follows. First, the 
quantum phase transition in iron-based superconductors 
with half-Dirac node described by half-Dirac fermion- 
Ising nematic field theory has stable fixed point A* = 
0, { v a/ {vf^f))* — with no relevant interaction away 
from it in the A— ^/(^f^f) parameter space and no non- 
trivial interacting fixed point. Second, on the dynamical 
properties of the quantum phase transition characterized 
in terms of response function describing quasiparticle- 
nematic coupling, we find that the system response ef- 
fectively has dynamical critical exponent 2 = 1. The 
nematic critical fluctuations lead to broadening of quasi- 
particle spectral function except in region near the node 
where we have 'Fermi arc'-like spectral ridge. Nematic 
fluctuations however do not affect much the degree of 
spectral peak anisotropy. However, the spectral ridge col- 
lapses in the vicinity of nematic QCP while it is well de- 
fined far away from QCP, both in ordered and disordered 
phases of nematic order. This signifies distinct quan- 
tum critical behavior from d-wave cuprates. This also 
emphasizes the distinction between the quantum critical 
(T = 0) behavior of quasiparticle Fermi arc and its clas- 
sical thermal phase transition counterpart, describable in 
terms of Kosterlitz-Thouless transition, with regards to 
the effect of order fluctuations (nematic Ising scalar order 
in the former and XY order in the latter). 
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We begin the formulation of the theory that give the 
above results in Section II by giving the low energy ef- 
fective action for this problem, considering the large- Nf 
version and computing the needed field theoretical quan- 
tities to proceed with the RG. Then we analyze the struc- 
ture of RG equations in Section III. We then discuss the 
nature of the quantum phase transition from response 
and spectral functions in Section IV. We end with dis- 
cussion of the results and connection with other works 
on the problem. 



II. LOW ENERGY EFFECTIVE FIELD 
THEORY OF SUPERCONDUCTING - ISING 
NEMATIC TRANSITION IN HALF-DIRAC 
NODAL IRON-BASED SUPERCONDUCTORS 

To describe half-Dirac nodes, we can start with 
standard Bogoliubov-de Gennes Hamiltonian describing 
Cooper paired electrons and approximate the Hamilto- 
nian around the nodes in momentum space [9]. We then 



construct the effective low energy field theory action for 
the superconducting state - Ising nematic order phase 
transition which consists of the action for half-Dirac 
nodal fermionic quasiparticles part, scalar Ising (ne- 
matic) order part and the fermionic quasiparticle-Ising 
order interaction part. The fcrmions describe fermionic 
quasiparticles 4*, * living in the vicinity of half-Dirac 
nodes in electron Fermi surface while the Ising nematic 
scalar field <f> describes the degree of lattice distortion 
from tetragonal lattice where (<fi) — to orthorhom- 
bic lattice where (<fi) ^ 0. The coupling between such 
nodal fcrmions with nematic order has been argued to 
be the most relevant coupling as the most effective pro- 
cess is where the nematic order field scatters the fermions 
around the nodes. 

In the actual physical situation, we have four fermion 
species for each of the pairs of (half-Dirac) nodes (11 and 
22), where we have spin up and down quasiparticles at 
one node and another set of spin up and down electrons at 
the partner node. In the technique, we generalize the 
spin up and down species into Nf "flavors" of fermions. 
The phenomenological field theory for this problem is 
then described by the following action 
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Sv= f 7^ V y]*i, (-iw m 7o + V/ft x 7l + ^^72)*l,a + *2,a(-*Wrr»70 +«/fc„7l + ¥^*x72)*2,a (1) 
J [Z7r > uj m a=l KF KF 



•* J ^ — i n „ — i 



n=l,2,a=l 



where 



*„, a (k,w„) = (c„ >0;+ (k,a;„),4 ia _(-k, -u n ), c n , 0i+ (k + Q, u n ), c\ a _ (-k - Q, -u n )) T (4) 



*n,o = *J,,«70 

and 

7o = n(3) T 3>7i = ^2 (^) t , 72 = ti(^)zt 2 (5) 

with index n = 1,2 representing the 11,22 pairs of 
nodes while a = 1,2, 3, .., Nf is the fermion flavor index. 
The momenta k, Q = kpx, k' = k + Q are as defined in 
Fig. 2. The r's are the usual 2x2 Pauli matrices. The 
details to obtain the fermion action are given in Appendix 
A. It is necessary to emphasize that Nf flavors do not 
imply spin Nf quantum number. The 4x 1 Nambu spinor 
defined in Eq. (4) contains the two Nf = 2 flavors from 
the two nodes in 11 or 22. 



Note that even though nematic order is normally rep- 
resented by symmetric traceless second rank tensor [1, 
7] due to the fact that inverting the "director" vector 
that characterizes the nematic state should give equiva- 
lent system, we write the Ginzburg-Landau action (2) for 
this theory in terms of Ising field because the free, non- 
interacting (without fermions) nematic phase transition 
falls within the Ising universality class. The coupling to 
fermions will in general able to change the universality 
class of the phase transition. 



Symmetry Consideration and Character of 
Nematic Order 



Nematic order is basically a Q = order that we have 
implicitly assumed to couple to spin singlet fermion bi- 
linear. It can be characterized in terms of particle-hole 
and particle-particle pairing correlators, 



(Ck'Ck') = -4k' 



(ck'C_k>) = (A^(-(cosk' x + 



• osk' ) + (l+cosk F )) + 1 



where we have represented our half-Dirac nodal 
anisotropic s-wave gap as A(k' x ,k' y ) = Ag(— (cosk' x + 
cosk' y ) + (l + coskF))- Here if is the overall superconduct- 
ing phase and we assume superconducting background 
through the entire order-disordered phases of the nematic 
order parameter. An ideal electron Fermi surface of iron- 
based superconductors with square lattice, without any 
symmetry breaking interactions, has C4(or C± v , which 
has extra symmetry of reflection with respect to ver- 
tical planes passing the central axis) symmetry. Note 
that both charge neutral(ones that do not depend on 
the sign of the gap) and charged (ones that do) observ- 
ables have this symmetry group, unlike in <i-wave super- 
conductors where only charge neutral observables have 
C^ v whereas charged ones have only Ci v . This symme- 
try group has four Id irreducible representations (s-wave 
with basis function f{k') = 1, d x i_ y i (cosk x — cosk' y ), d xy 
{sink' x sink' y ) , and g {sink' x sink' y {cosk' x — cosk' y )) and one 
2d irreducible representations p(sink' x , sink' y )[\\]. With 
such pairing gap, it is clear that A k > = 0, = <fi 
(where <fi is real field) breaks C^ v symmetry to Ci v . This 
is nematic order. So, symmetry arguments that char- 
acterize how nematic order couples to our half-Dirac 
nodal superconducting state suggest that nematic or- 
der parameter has s-wave symmetry where time-reversal 
symmetry is unbroken but the point group symmetry 
is reduced from Ci v to Civ with basis function sim- 
ply a unity function. The nematic order can be po- 
larized along x or y direction. This extra s-wave com- 
ponent cf> however will eliminate the half-Dirac nodes 
for <j) > or change each of the half-Dirac nodes to 
two full Dirac nodes for (f> < 0. We can however have 
order parameter which breaks C± v to Cq, v while keep- 
ing the half-Dirac nodes intact which clearly must be 
higher order modes. From the 5 irreducible representa- 
tions of Ci v point group symmetry, d xy wave with ba- 
sis function sink' sink',, can deliver such half-Dirac nodc- 
preserving C^-Xo-C^v symmetry breaking order param- 
eter. To be precise, we can add on top of our anisotropic 
gap A(&4, k' y ) = Ag(l + coskp — {cosk' x + cosk' y )) an or- 
der parameter B^ = <psin(k' x — k c )sin(k' y — k c ) where 
< k c < kp/y/2. The resulting gap for <f) < is shown 
in Fig. 3 while for <p > the gap is Fig. 3 rotated by | . 
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FIG. 3 Cav to Civ symmetric breaking of electron Fermi 
pocket gap with half-Dirac nodes. 

The structural tetragonal to orthorhombic phase tran- 
sition can easily be seen by rotating the gap function by 
? to align with the lattice x and y axes. 

The order parameter fluctuations when coupled to 
fermions will be relevant only when the wave vector Q 
carried by the order parameter is also the wave vector Q 
that connects the two nodes between which the fermions 
are scattered, by momentum conservation[ll]. That is, 
the order parameter will scatter the fermions effectively 
and efficiently only when the momentum it carries is 
transferred entirely to the fermions. Otherwise the scat- 
tering is a virtual process which merely renormalizes the 
coupling constant without making a fundamental change 
in low energy theory. In this case, nematic order corre- 
sponds to Q = which means our theory consider scat- 
tering of fermions living in the vicinity of the same half- 
Dirac node by the nematic order fluctuations. Scattering 
by Q = {2kp, 0) or Q = (0, 2kp) should be described by 
some type of density wave but if we want to couple this 
density wave with fermions, the theory should describe 
coexisting superconducting and density wave phases. Ex- 
perimentally as shown in Fig. 3 of Ref. 16, the super- 
conducting phase does not coexist with the density wave 
phase but they do border each other, even though the 
Ising nematic order, where we have orthorhombic crys- 
tal structure, penetrates deep inside the superconducting 
regime, as illustrated earlier in Fig. lb). We therefore 
essentially consider Ising-nematic quantum phase transi- 
tion within the background superconductivity. 



B. The Effective Action of Ising Nematic Order 
Parameter 



In the large-TV/ expansion, with the physical case 
Nf = 2 corresponds to spins up and down, we gener- 
alize the two spin polarizations of S = \ fermions into 
Nf "flavors" of fermions. We compute the effective ac- 
tion for bosonic nematic order parameter field S e f / [<p\ by 
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formally integrating out the fermion tp, ip from the orig- 
inal full action S[ip,</>]. This will give rise to nonlocal 



logarithmic term containing the fermion-nematic order 
field coupling constant A and can thus be expanded per- 
turbatively in powers of A. 



S eff [4>] = S+-Nf ]T In(Det[- 



+ lt n F .k 7l + .k) 2 72 + \<hi]) 

vakf 



(7) 



Rescaling <j> $ and r — > rNf\ 2 and retaining the 
surviving terms, we have 

S eff [4>]=N f Sl ff 

where 



Sl f/ = i <?x\r(t> 2 - V In(Det[-iu ml0 + ^£.k 7l + -^(^« .k) 2 72 + <fryi]) (8) 
J 2 J-f v A k F 



n=l,2 



where the (d T cj)) 2 , (V0) 2 , u<fi 4 terms vanish as we take 
the Nf — > oo limit. We will only consider the quadratic 
correction terms in <p while assuming that the renor- 
malized quartic terms will remain sufficient to stabilize 
the effective Ginzburg-Landau type effective theory for 
(j> without explicitly considering its renormalization. The 
effective action to quadratic order for (j> can be written 
as 



d 3 k 



where 



r 2 (fe) = U 2 (k x ,k y ,uj) +U 2 (k y ,k x ,uj) 



(9) 



(10) 



The polarization function Tl2{k x , k y , uj) is given by the 
Feynman diagram in Fig. 4 



k+p 




FIG. 4 The Feynman diagram for polarization function 
which contributes correction to quadratic part of <j> effective 
action. 

Using the representation (4) for Dirac 7 matrices, the 
expression for the polarization function which is the cor- 
rection to the boson propagator due to fermion is given 
by 



Tl2{k x ,k y ,uj) 



with 



d 3 p 



TrfriGvip^Gvip + k)} (11) 



G 1 (p) = -lUJ-f + VfPx-fi + ^^P 2 y l2 (12) 

where £ = and we have focus only on node 1 
as example. Despite its lengthy form, this polarization 
function is still even under time reversal of external mo- 
menta (k x ,k y ,uj) —} (—k x ,—ky,—ui) which can be veri- 
fied by direct inspection. Also, by power counting, this 
expression must have dimension one in external momenta 
[Il2(fc)] = k 1 but is an even function of k = (k x , k y ,oj). 




FIG. 5 The profile of polarization function Il2(k x , k y , u) = 
— O.OOlAo;) computed numerically with vf — va = 0.067eF^ 
where a = 0.1 is the lattice spacing, with one unit of length 
= 100^4°. One unit of momentum Ak = ^ and frequency 
Auj — 0.271meV . Here we chop off the tip to show the sym- 
metry of the profile. The function has even parity in (k, w) 
and C4 symmetry. 
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We compute numerically the value of Il2(fc) in Eq.(8) 
as function of the two external momenta at a fixed ex- 
ternal frequency and the result is shown in Fig. 5, which 
gives us some rough idea of how this quantity varies as 
function of momenta. The most important observation 
from Fig. 5 is that the polarization function profile has 
C 4 symmetry and that it is peaked around k' = (0, 0) 
with respect to the center of electron pocket. 

We also compute the fermion self-energy and vertex 
corrections to nematic scalar field action needed for the 
following sections, the details of which are given in Ap- 
pendix B. Regarding the nature of quantum phase tran- 
sition, we will show that the quantum phase transition 
between Ising nematic ordered and disordered phases, in 
the presence of half-Dirac fermions in the background 
superconducting state, is second order phase transition 
with the non-interacting fixed point of Dirac fermions. 

III. RENORMALIZATION GROUP 
EQUATIONS STRUCTURE OF THE THEORY 

One thing that can be noticed from fermion action (1) 
is that the momentum is linear in k x but quadratic in 
k y . If we want this action to be invariant under rescaling 
of momenta and frequency, the scaling dimension for k x 
must be twice of that of k x . But here for the purpose 
of RG calculations we choose to use the same scaling 
dimension for k x and k y . Further, rather than using a 
dynamical critical exponent z for the scaling of frequency, 
we use the same scaling dimension for u> as that for k x 
and k y . The reason is RG will naturally fix the scaling 
dimension for to implicitly inside the coefficients of RG 
equations. Therefore, we write the scaling of momenta, 
frequency as follows; 



here with Yukawa type of coupling, it can be easily 
checked that the scaling dimension of Yukawa coupling 
constant is dim[X] = \ which shows that the coupling to 
nematic order parameter is a relevant perturbation to the 
( fcpt F )* = ^ nxe d point associated with the half-Dirac 
nodal quasiparticlcs of the iron-based superconductors. 
The mass r term on the other hand has scaling dimen- 
sion dim[r] = \ = 2 at tree level which shows it is also 
relevant parameter and must thus be evaluated in the RG 
analysis, the RG flow equation of which can in general 
be written as 

d i = lr = (2- Vb )r (15) 

where the generally nonvanishing anomalous dimen- 
sion r\ h of bosonic scalar field <j) enters. As a comparison, 
one may prefer to use a rescaling of momenta and fre- 
quency that leaves the fermion action invariant, 

k x = k x e 



u = uj'e- 1 (16) 

Again considering renormalization of the time deriva- 
tive term in and respectively gives scaling dimen- 
sions dim[$(k)] — — | and dim\$>{k)] = — \. This implies 
field scaling of the form 

$(fc) = $'(A0e-/o 1) 



k x — k x e 
ky = k y e 

lu = uj'e- 1 (13) 

The scaling of bosonic and fcrmionic scalar fields is ob- 
tained by considering renormalization of the time deriva- 
tive term in and S^p respectively which gives scaling 
dimensions dim[<&(k)\ — — | and dim[^(k)} = —2. This 
implies field scaling of the form 

$(fc x , ky,u) = &{k x ,k v ,uj)e- Si 

*(k x , ky, u) = *'(fc x , k y , u)e~ So <«'(2-¥) (14) 

By similar simple power counting at tree level, where 
interaction is assumed to be nonsingular, as is the case 



yf(k) = ^i'{k)e-So dl '(l-^) (17) 

We still arrive at the same conclusion that the 
Yukawa coupling is relevant perturbation away from 
the ( k py F )* = fixed point with scaling dimension of 
dim [A] = |. 

We investigate the structure of the RG equations 
by starting from the logarithm of cut-off derivative of 
fermion self energy and vertex functions rather than us- 
ing the usual momentum shell transformation. The RG 
equation should be of the form given by 

A JJ = (Co 7o+C 1 7i)^+(C 1 7o+C 1 1 7i)fcx+C 2 (fc,) (18) 



A£=C 311 (19) 

The k y part is written to indicate that it has some un- 
known nonlinear dependence on k y , unlike the k x and lu 
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parts. We also note that the inverse propagator is not 
homogeneous function of u>, k x and k y anymore and this 
has tremendous effect on the complexity of the analy- 
sis. Before proceeding further, we first find the relations 
between the C's coefficients and the anomalous dimen- 
sions rjf and r/i,. One caveat before doing this is that the 
RG equations (14) and (15) are 4x4 matrix equations 
rather than scalar equations. The relation involving rjf 
is obtained by fixing the kinetic time derivative part in 
while that for rjb is obtained by fixing the Yukawa 
coupling term. The latter is done because the nematic 
order action does not have kinetic term upon inte- 
grating out the fermions, as can be seen in Eq. (6) which 
has a non-local term given by the logarithmic part. This 



non-local term will overpower the local terms of S^, ex- 
cept the r(fi 2 term. This can be seen by considering the 
scaling dimensions of each of the terms in S^. Earlier 
we obtained dim[$(k)] = — | or dim[<j)(x, y, t)] — | and 
this implies dim[r] = 2 and dim[u] = 1. We see that 
the quartic coupling u is less relevant than the quadratic 
coupling r and so the former can be omitted. The deriva- 
tive terms on the other hand will also have diminishing 
effect because they are marginal terms. We are therefore 
allowed to retain only the rcf> 2 and omit the other terms 
in S^. 

To find the relations between the C's coefficients and 
the anomalous dimensions rjf and rjb, fixing the kinetic 
time derivative part in S\$, we have 



Vf = 7o~ 1 Co 7o + C 1 7l 



C\ - C\ 1 

o c\ + cl 
o c\+cl 
ooo c\-cl 



(20) 



r 



which demonstrates a "pair- wise" pattern associating 
fermions living at the same nodes (fermions 1 and 4 at 
node 1 (or 2) and fermions 2 and 3 at node l(or 2). Wc 
will also see this pattern later in the structure of RG flow 
equations for v F and da which has crucial significance as 
hinted in the abstract. If we argue that the anomalous di- 
mension rjf is a universal characteristic of fermion which 
should be independent of the where the fermion lives in 
momentum space, then wc have to conclude 



CO /^fO /"rO 
— Ls 1 — O 



C? = Vf 



and therefore 



Cg = »7/,C?=0 



(21) 



Following similar procedure to rjb by fixing the Yukawa 
coupling, we have 



rj b = 1 + 26*3-277/ 



(22) 



Here, C,t)/<1 but rjb <~ 1. 

The observation that we only have one equation for 
boson field anomalous dimension fits perfectly with the 
fact that we only have one real scalar boson field. We 
derive these RG equations for coefficients C's by evalu- 
ating directly from the expression for the self en- 
ergy T,(k,u>). Unfortunately, since the propagator is not 
homogeneous function of momenta and frequency, the 
resulting equations will depend explicitly on the cut-off 
momentum A. We can use hard cut-off since this is the 
simplest possible cut-off function and hope this choice 
will minimize the cut-off dependence of the final results. 



Computation of C's as functions of v F and ^a can only 
be done numerically. Rather than pursuing the detailed 
numerical evaluation of the flow, we instead analyze the 
structure of RG flow equations of the theory and make 
qualitative but yet definitive conclusions about the crit- 
ical properties of the system. 

The RG equations for Fermi and gap velocities are de- 
rived by considering the renormalization of terms that 
contain these velocities in the fermion action . Fixing 
the VFk x term, we have 



dVF 

~df 



(-Vf - Chi So + C\)v F 



(23) 



where v F = (v F1 ,v F2 ,v F3 ,v F4: ) T 
whereas from renormalization of 



8-uaP 



- term we have 



d,VA 



= (~Vf ~ (G 2 7 2 " 1 7o + C 2 ^i)>a 



(24) 



where wa = (^Aij v A2, v A3, v aa) T ■ We only consider 
the leading linear order terms for the RG equations here. 
Using the 7 matrices as defined in (4), we have 



/ dvpi \ 
dl \ 



1 dv F4 , 

\ dl ' 



dvf 

dl 
dv F2 
dl 

dVF3 

dl 

dVf4 
ST 





M a 





1 








M 
















M 
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M a 





fv F1 \ 


1 


V F2 




1 VF3 




\ VF4 J 



(25) 



where 



M a = - Vf - C\ + C°, = -rjf - C{ - C? 
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dl } 

a«A2 






dl 


"I 










dl 




V 


*UA4 
dl / 






We have RG flow equations for Vp and i>a involving co- 
efficients C's that are themselves functions of vp and i>a- 
We can see here that each of the four fermion components 
of the 4x f Nambu spinor has its own RG equations for vp 
and t>A- This means that in general vp and v& will flow 
differently between these four fermion components under 
RG. This feature arises purely from the symmetry prop- 
erty of anisotropic gap of electron pocket which is clearly 
different from c?-wave cuprates case where one RG equa- 
tion for vf and another RG equation for va are sufficient. 
If we consider the hypothesized fixed point ( v t )* = 0, 
which necessarily requires («a)* = but (vp)* ^ 0, from 

and C\ = —rjf at the fixed 



^f=0we will have Cft 



point which set the final condition for the RG flow equa- 
tion for the coefficients Cq as Cq(v a = 0,v F ^ 0) = 
and C\ as C\(v* A = 0,v* F ^ 0) = —r\f- This fixes the 
final condition for both coefficients in the RG flow for 
vp. Considering the = at fixed point, we obtain 
no new information but the confirmation that the RG 
flow equation is consistent at the fixed point. This 
consistency at least offers some posteriori justification of 
our fixed point assumption. 

If we denote the four fermions constituting the 4x1 
Nambu spinor Eq. (4) by A, B 7 C and D respectively, 
as depicted in Fig. 2 so that \t„ ja = (A, B,C, D) T , 
looking at the structure of the flow equations given in 
equations (22) and (23), we can see that the flow equa- 
tions for fermion A have equivalent structure to those 
of fermion D while those of fermion B are equivalent in 
structure to those of fermion C. This means that the 
RG flow associated with the spin up nodal fermion at 
(kp 7 0) is the same as that of spin down nodal fermion 
at (A)f,0). Likewise, the RG flow of the spin up nodal 
fermion at (— kp,0) is identical to that of spin down 
nodal fermion at (— hp, 0). This agrees with the expecta- 
tion that the Fermi velocity of two fermions living at the 
same half-Dirac node must flow identically but this also 
shows a new important observation that the Fermi veloc- 
ity at the two half-Dirac nodes can in general flow differ- 
ently under RG. One direct physical implication of Fermi 
velocity flowing to different values between fermions at 
different nodes is that, since e k = v F k x + 0(k 2 ) and 
v A = V fe A(fc)| k=(feF;0) = the originally C 4 symmet- 
ric electron pocket gap is deformed under RG flows (where 
nematic order comes into play implicitly via its coupling 
to nodal Fermions) and thus the C 4 is broken(Fig. 3)). 

The relevance of nematic-quasiparticle Yukawa cou- 
pling at tree level led to the suggestion of the existence 
of a non-trivial fixed point A* ^ along the coupling 
A axis [2]. This fixed point thus has (^)* = along 
the other axis. We have shown earlier that with half- 



Dirac spectrum, nematic order is still relevant coupling 
despite with slightly different scaling dimension. There- 
fore at tree level one would naively expect that A will 
flow away from the noninteracting Gaussian fixed point 
A* = 0. However, if we consider one loop diagrams for 
renormalization within minimal subtraction scheme and 
dimensional regularization, simple power counting shows 
that the all the diagrams have no divergence. This shows 
that the Yukawa coupling has vanishing ji function and 
is thus irrelevant coupling at one loop level, despite being 
relevant at tree level. 

One of the signatures of the proposed Ising nematic 
phase is its effect on the superconductor quasiparticle 
spectral function. In particular, quantum critical ne- 
matic fluctuations are expected to damp the quasipar- 
ticle spectral function, resulting in significant change in 
the spectra. In d-wave superconductor, it was shown 
[2] that quantum critical nematic fluctuations overdamps 
the quasiparticles in direction normal to the Fermi sur- 
face while weakly damps them in direction tangential to 
the Fermi surface around the node, resulting in highly 
anisotropic spectral function in the form of very narrow 
wedge. We will show soon that in the case of iron-based 
superconductors with anisotropic gap in electron pocket 
the situation is rather different. 



IV. THE NATURE OF THE QUANTUM PHASE 
TRANSITION FROM DYNAMICAL PROPERTY 
OF THE THEORY 



In this section, we will investigate the nature of this 
quantum phase transition [18,19,20] from the system's 
critical dynamics. The interaction between half-Dirac 
nodal quasiparticle and nematic Ising order field can be 
interpreted as some kind of "damping" process where 
Ising order parameter <f) decays into half-Dirac nodal 
quasiparticles One can then write a phenomenological 
field theory characterized by a dynamical critical expo- 
nent z that measures the degree of the damping of the 
half-Dirac nodal quasiparticles by the nematic order pa- 
rameter fluctuations. To quantify this dynamical process, 
we consider the finite temperature version of polarization 
function T2(k, a;„) to Ising nematic scalar ordering field 
action due to half-Dirac nodal fermions (8). Since polar- 
ization function can also be interpreted as susceptibility 
Xo(k,uj n ) describing the system response, we can also 
use it to understand the nature of quantum critical phe- 
nomena[27,28] from this dynamical property around the 
quantum critical point at optimal doping level x = x Cl 
assuming that the half-Dirac node is intact by control 
of appropriate parameters. The susceptibility function 
Xo{k,uj n ) applicable to the finite temperature crossover 
around QCP is given by 
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(27) 

where according to Eq. (8) we have to include the con- 
tribution of both pairs 11 and 22 of nodes. We are not 
particularly interested in the explicit form of the above 
expression but rather, in its form upon lowest order ex- 
pansion in powers of k x , k y and uj n where k x , k y are mea- 
sured from the node as shown in Fig. 2. It can be shown 
that 



Xo(k x ,ky,u n ) = xo{0,0,0) + aul + b(kl + kl) + .... (28) 

which shows the presence of kinetic term w 2 ver- 
sus kl,ky and thus implies dynamical critical exponent 
z = 1. We can contrast this with the Fermi liquid-spin 
density wave transition[18] for example which has low 
energy susceptibility of the form xo(k,u n ) = Xo(0, 0) — 
Ci|w„| — c 2 k 2 + .... where the \co n \ indicates the damping 
of order parameter fluctuations due to due to the cou- 
pling to nodal fermions at the nodes connected by the 
spin density wave ordering wavevector. Note that we ob- 
tain the above result Eq. (34) even without taking the 
asymptotic limit va — > or making a priori assumption 
that we have fixed point at (ua/(i>f&f))* = 0. This 
therefore independently agrees with our result that the 
field theory Eqs. (2,3,4) has fixed point at (j^)* = 
because then the fermions behave like those described by 
action ~ ^^ a {-iuj m T z + VfKt*)^!^ + ^2, a (-ioJ m T z + 
VFk y T t )^2,a asymptotically in the limit v& — >• which 
has z = l(In calculations involving fermion internal mo- 
menta we cannot take «a = though) . 

We observe here that despite having Lorentz-symmetry 
breaking anisotropic dispersion, in the low energy limit 
our Ising nematic field theory behaves as undamped sys- 
tem characterized by z = 1. This suggests the Ising ne- 
matic field fluctuations are effectively undamped by its 
coupling to half-Dirac nodal fermions. This is only pos- 
sible if (or rather, this implies that), at nemating quan- 



tum critical point, the quasiparticle density of states is 
very low or actually vanishes. We will shortly show more 
quantitatively that this is indeed the case, by directly 
computing the quasiparticle spectral function. Nematic 
critical fluctuations however do not necessarily change 
much the degree of anisotropy of quasiparticle spectrum, 
but rather, the quasiparticle spectral weight itself. This 
is to be compared with the result for the d-wave case [2] 
which found a very strong damping of nematic order pa- 
rameter by quasiparticle excitations in such a way that 
the spectral function of the quasiparticle is significantly 
broadened everywhere in BZ except at narrow 'wedges' 
around Dirac nodes where the spectral function acquires 
a very anisotropic Fermi arc shape, being very narrow 
in direction perpendicular to the Fermi surface and very 
long in direction tangential to it. In other words, ne- 
matic critical fluctuations strongly enhance the velocity 
anisotropy of the <i-wave nodal fermions. This later con- 
clusion is however, as pointed out by Pclissetto et al.[21], 
very dependent on the use of nematic scalar field and tree 
level power counting result which casts the Yukawa cou- 
pling between nematic order and nodal quasiparticle as 
relevant. Higher loops and field theorist renormalization 
calculation of scaling dimension of the coupling constant 
of composite spin density wave rather than Ising nematic 
to nodal quasiparticles puts the coupling as marginal 
rather than relevant. 

We now compute the half-Dirac quasiparticle spectral 
function to see more quantitatively the effect of nematic 
critical fluctuations. The spectral function is given by 

A(k, u) = -2sgn{u)Im[G ab (k, «)] (29) 

where G a {,(k, to) is the nonvanishing (a, b) — th element 
of the renormalized single quasiparticle Green's function 
given by G _1 (k,w) = C 1 (k,uj) — E(k,w). We consider 
the 4x4 matrix version of our theory with Go given in 
Eq. (10) and the self energy E(k, u>) can be decomposed 
as E(k,w) = E a 7o + E b 7l + E c 72 , the details of which 
are given in Appendix C. We compute the polarization 
function numerically and hence the self energy as well 
and input the result to the following expression for the 
spectral function, 



w + 7ro(E a (k,w)) 
(Jro(£«(k, w)) + w) 2 + (Z\k,co)-v F k x y + (£<=(k,w) - ^klf 



We present here the dependence of quasiparticle spec- tral function on velocity ratio v^/vp shown in Fig. 6. 
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FIG. 6 The quasiparticle spectral function A(k, uj = 
— O.OOIAoj) centered at the node k' = (k,F,Q) with respect 
to the center of electron pocket, as function of anisotropy 
vf/va- We use as parameter values; one unit of momentum 
Ak — j^, the lattice spacing a = 0.1, with one unit of length 
= 100A° and frequency Au = Q.271meV and at vf/va = 1, 
we use vf — va = 0.067eV-. 



We note that in the limit of vp ^> va, we have ex- 
tremely anisotropic spectral weight along one of the two 
(e.g. x) orthogonal axes. As the ratio between the two 
velocities is tuned to order of unity, the anisotropy de- 
creases but we still have relatively anisotropic spectral 
function. At vf ~ "a the spectral peak is actually still 
very anisotropic with 'Fermi arc'-like shape rather than 
round one. It takes some large ratio vf/va > 1 in or- 
der to achieve a roughly round spectral peak. At the 
other limit of the anisotropy where va 3> vp, we have 
extremely anisotropic ridge perpendicular to the ridge in 
the opposite limit. This behavior clearly derives from the 
non-Lorentz symmetric form of fermion action. 

The extent of the effect of nematic critical fluctua- 
tions on quasiparticles is best investigated by consider- 
ing the dependence of spectral function on the strength 
of quasiparticle-nematic Yukawa coupling A. We expect 
that, in terms of anisotropy, the coupling to nematic or- 
der does not change much the anisotropy of the quasi- 
particle spectral function. The effect of nematic order 
is however more towards the magnitude of spectral peak 
and we expect that there could be nontrivial change in 
the profile of spectral function near some critical value 
A = A c with regard to the presence or absence of Fermi 
arc. This can readily be expected from the fact that 
the effective action S\^\(j)\ in Eq. (8) can eventually be 
written as, 



d 3 k 

(2tt) 3 



(r-r(A)Mfc)| s 



d 3 k 
(2tt) 3 



(A 2 -A 2 )|<Kfc)| : 



(31) 



where the mass parameter r ~ A 2 corresponds to 
the critical coupling while r(A) ~ A 2 comes from the 
the Yukawa nematic-quasiparticle coupling. Intuitively 
speaking, for A -C A c (deep inside disordered state of <j>), 
we have (4>) — and so that in computing the renor- 
malized fermionic quasiparticle propagator we can sim- 
ply substitute ((f)) = for <f>. Likewise, for A A c (deep 



inside ordered state), we have (4>) ^ and we can sim- 
ply substitute this (<p) for 4>. In both cases, the quasi- 
particle spectral ridge is well defined. Things are less 
clear in intermediate region near A c as to what happens 
to the spectral peak. To investigate the physics in this 
regime and also to verify the simple intuitive picture 
in deep ordered and disordered phases, we again com- 
pute numerically the spectral function but this time we 
vary the Yukawa nematic-fermion coupling strength A. 
We find that spectral function has sharp ridge in the (f>- 
disordered state (A < A c ) and much heavily broadened 
spectral peak in ordered state (corresponding to mas- 
sive fermions) but the spectral peak collapses in the in- 
termediate region, which directly corresponds to critical 
point A c in S 1 ^^ [</>]. This can actually be checked ana- 
lytically. From Eqs. (37) and (C2) in Appendix C, it is 
clear that j4(k, ui) vanishes at a critical value A c , which 
in the limit of r 3> A 2 , is approximately given by 



r M 



d 3 p 



(32) 



(2^ [( W + 0)2 +t ,J,p2 + ^2 p 4] 



with j] = and we have focused on the nodal point 
(k x = k y = 0) at which the peak is centered. This re- 
sult clearly quantitatively confirms our earlier intuitive 
prediction. This T = critical point shows different be- 
havior with regard to spectral peak than what occurs 
in thermal phase transition in cuprates where Fermi arc 
emerges as the temperature (here the Yukawa coupling 
A acts as temperature with A c — A ~ T — Tkt) is in- 
creased above Tkt i n a. Kosterlitz-Thouless type of tran- 
sition when analyzed with XY model [29]. In the QPT 
described by Eq. (37), in both deep ordered and dis- 
ordered regions of <f> field, the fermionic quasiparticles 
behave as free fermions (but being massive and massless 
respectively) and are well defined in these two. Only in 
the vicinity of QCP the nematic quantum fluctuations 
are able to disrupt the definition of quasiparticles. This 
underlines the difference of the physics of quantum and 
thermal phase transition in these high T c superconduc- 
tors and at the same time provide validity check to our 
theory. 
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V. DISCUSSION AND SUMMARY 

We have emphasized the analysis of RG structure of 
the low energy effective theory of half-Dirac iron-based 
superconductors to deduce its nematic quantum critical 
properties. Simple power counting at tree level shows 
that nematic order is relevant perturbation away from 
the ( k v * )* = fixed point of the half-Dirac nodal quasi- 
particfes. However, all one-loop diagrams in the context 
of dimensional regularization plus minimal subtraction 
scheme of field theorist rcnormalization have no diver- 
gence which suggests that at one loop level the nematic 
order is irrelevant coupling and the theory thus has non- 
interacting fixed point and that { k v J^ ) is irrelevant pa- 
rameter. The RG flow diagram therefore has stable flow 
lines both along the Yukawa coupling A axis and as well as 
the ( ) axis. The theory therefore does not have finite 
coupling fixed point which is the case for d-wave cuprate 
superconductor. The C4 symmetry breaking leads to pe- 
culiar nature of relevance(or irrelevance)of the RG pa- 
rameters around the electron pocket where vf is found 
to be marginal coefficient while , " A can be marginal or 
irrelevant depending on which of the two equally valid op- 
tions of momenta-frequency rescaling scheme we choose. 
The k x O k y equivalence of C4 symmetry is broken by 
this anisotropy and this suggests that the gap deforma- 
tion instability associated with nematic order is the phe- 
nomenological physical picture of structural phase transi- 
tion in half Dirac nodal iron-based superconductors. The 
whole analysis assumed that the half-Dirac nodes remain 
intact all along during RG flow. The possibility for this 
to be the real situation is strongly supported by recent 
work in Rcf. 26 where it was shown the such kind of 
nodes is guaranteed to exist as long as the strength A 
of hybridization between the two electron pockets (the 
interpocket hopping term with momentum (it, it, tt)) is 
less than some critical value A c . Our theory is therefore 
valid within finite regime in parameter space rather than 
only at a critical point that can only be achieved by fine 
tuning. 

In describing the physical mechanism of structural 
phase transition in high T c superconductors, the nematic 
phase couples most relevantly to the quasiparticles while 
at the same time couples to lattice distortion which mea- 
sures the degree of structur deformation. The nematic 
phase can be electronic (charge) or spin nematic phase 
and is still currently actively investigated experimen- 
tally [24]. The gap anisotropy in iron-based superconduc- 
tors is believed to be determined by orbital content [25] 
of the electron Fermi pocket and this directly suggests 
connection of structural phase transition to orbital or- 
dering as orbital ordering is driven by redistribution of 
occupation density of d xz and d yz orbitals in (tt, 0) and 
(0, tt) electron pockcts(in extended Brillouin zone). We 
see therefore a self-consistency of picture of the physics 
of structural phase transition as related to the symme- 
try breaking of anisotropic gap on the electron pocket 
with the proposal of orbital ordcring-driven structural 



phase transition. It was also shown that in orbital order- 
ing picture [17] the structural phase transition is in the 
Ising universality class and can be described by effective 
Hamiltonian 



Hspt = -Jspt M - M i ( 33 ) 
{i,o) 

which is not inconsistent with our Ising nematic order- 
ing picture and idea that nematic ordering can be related 
to orbital ordering. 

Another important support for nematic-driven struc- 
tural phase transition[22] is provided by the order of 
phase transition. In the 122-family, the structural and 
magnetic ordering phase transitions become a single first 
order phase transition while in 1111-family, they become 
two separate second order phase transitions. Eventually 
the idea of nematic phase offers very good explanation 
of this kind of behavior [23]. In magnetic ordering, the 
order parameter is characterized by vector while in ne- 
matic phase the order parameter is director, which is 
headless vector which we have represented with Ising ne- 
matic scalar field because the phase transition is in Ising 
universality class. While the magnetic ordering phase 
transition is second order, the nematic phase transition 
is known to be very weakly first order. This suggests that 
in 122 compound, the nematic ordering is very dominant 
and controls the magnetic ordering so that they become 
a single first order phase transition. In 1111 compound 
on the other hand, the nematic and magnetic orderings 
are equally strong and occur in sequence where the cubic 
term in the Ginzburg-Landau free energy description of 
nematic order diminishes and we thus obtain our effec- 
tive field theory description Eq. (15) and second order 
phase transition for both ordering. 

The very existence of nematic phase in high T c super- 
conductors including iron-based superconductors is ac- 
tually still a very much debated issue. Several impor- 
tant fundamental questions include, if such phase does 
indeed exist, what kind of nematic phase it is precisely, 
electron nematic or spin nematic? It is then interesting 
to ask the differences in physical consequences of these 
two (or other) possible forms of nematic phase. Also, 
if indeed this phase exists, where does it exists in the 
phase diagram? In the paramagnetic tetragonal phase at 
T > T s or in the magnetically ordered phase T < T N 
or in between the two or instead it spreads over different 
regimes? Latest experiment [24] suggested the presence 
of an intrinsic spin nematic fluid even in the isotropic 
tetragonal phase T > Ts of BaFe^As^ based on the ob- 
served anisotropy of spin wave excitations that persist 
well inside this region in the phase diagram. It is an 
open problem left for future work to construct a theory 
based on spin nematic fluid that can explain structural 
phase transition in iron-based superconductors. 
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nodes respectively. The kinetic energy is measured rela- 
tive to Fermi energy ;(k = ek — Sf where k = K + q with 
2K = Q as depicted in Fig.Al). Simple inspection on the 
geometry of the gap suggests that we have £k = — Ck an d 
A(k - Qi) = A(k ) = A(-k). We assume that the the- 
ory has standard BCS phenomenology with anisotropic 
s-wave symmetry at the electron Fermi surface with well 
defined, long-lived fermionic quasiparticles. We aim to 
write the fermion action as 



Appendix A: Details on the Fermion Action 

We will give the details here the construction of 
fermionic quasiparticles action. We consider electron 
pocket with anisotropic gap associated with the fermionic 
quasiparticles, as illustrated in Fig. Al). 



d 2 k 



S = I 72^ T E *„, (k,o; n )M*„ )0 (k ) a; n ) 

J \ J Un n=l,2,a=l 

(Al) 

where ^(k,^) is as given in Eq.(4). The general 
quasiparticle action at T ^ takes the form 




FIG. Al)Electron Fermi surface with critical half-Dirac 
node 

The 'nesting' wave vector in this case is Q — (±2fcp, 0) 
or Q — (0,±2fcp), corresponding to 11 and 22 pairs of 



d 2 k 

S= I , n ^ T} [(iuj n -( k )cl r (k,uj n )c cr ('k,u} n ) 



(2tt; 



2 i 

<T,LU n 



- - A fe 4(k, w„)4(-k, -w n ) + h.c. + 0(c 4 ) (A2) 

We thus just have to find the elements of matrix M 
by matching the corresponding terms in the two forms of 
action. Before doing that, we need to separate the sum 
over momenta into those over k and those over k — Q. 
Doing this, it can be checked that the action takes the 
form 



d 2 k v 



N f 



*„,„(k,w n )(io; n - CfeMi - A fc M 2 )*„ !Q (k,a;„ 



(A3) 



n=l,2,a=l 



where 



Mi = 



( i 

-i 

i 

Vo -i 



10 
10 
1 
10 



(A4) 



as example. The resulting effective action is not Lorentz 
invariant even if we try to rescale the coefficients in front 
of the operators in the bracketed terms. The first two 
terms in the bracket can however be treated as if they 
form a half-Dirac action and we can thus write the Dirac 
representation for these two terms which demands us to 
construct two anti-commutating 4x4 Dirac 7 matrices. 
We can of course define another 7 matrix, call it 72, for 
the third term but that is not bound to satisfying anti- 
commutation relation with the first two 7 matrices. We 
choose the following representation; 



Note that the form of M\ and M 2 is fixed by the form 
of kinetic energy and the gap symmetry A& respec- 
tively, and is therefore unique. These are expanded to 

SAk 2 

lowest order as ~ vpk x and A& ~ , y for node 1 



7o 



10 
-1 
10 
0-100 
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7i 



72 



/ 10 
1 
-10 

V o -loo. 



/ 1 \ 

0-10 

10 

V -l o o o / 



(A5) 



The choice of 71 for the "Yukawa coupling" for inter- 
action action in equation 3 is derived based on symmetry 
consideration as follows. The fcrmion-boson coupling ac- 
tion must clearly preserve C4 symmetry of the theory 
since the theory must be invariant if we rotate the mo- 
mentum space by Under this rotation, the operators 
in the 4x1 Nambu spinor in Eq.(4) transforms as 



Cl,a+(k,U>„) C 2 , a +(k, U n ) -> Ci, a _(-k, -U n ) 



c l,a-(- k .- W «) c 2,a-(- k '- W n) ~> C f la+ (k,u;„) 



ci, 0+ (k- Qi,w„) -)■ c 2 , a +(k- Q 2 ,^n) Ci, _(-k+ Qi,-w„) 



4,a-(- k +Ql.- W n) -> 4,a-(- k + Q2> ~> 4,a+( k ~ Ql^„) 



(A6) 



With lengthy but otherwise straightforward algebra it therefore we adopt this particular choice of Yukawa 



can be verified directly that 



coupling. 



N f 

-SW = A / ^ / dT J2 0*n,a7l*n,a (A7) 
^ n=l,o=l,2 

is invariant under this C4 transforma- 
tion, while other choices such as S^^ = 

\Jd 2 xJdr J2nii,a=i 4>^n,alo^n,a are not. Invari- 
ance under C 4 symmetry is the minimum symmetry 
requirement that must be satisfied by the action and 



Appendix B: Self-energy Correction to Fermion 
Propagator and Yukawa Vertex Correction 

We consider a theory of fermions with Nf flavors, writ- 
ten as 2 x 1 spinor which means we consider only one of 
the two nodes in each of 11 and 22 pairs, interacting with 
real scalar field through Yukawa coupling. 



Nf 

S *= I T^^EE^H^* +v f k x T t + ^kliTy)^ ha +^ 2 , a (-iiJ m r z +v f k y r t + ^k 2 jTy)^ a (Bl) 

J ^ Zn > u> m a=l kF KF 



and = ^J, t*, where oj n is the Matsubara (imag- 

/l 2 1 2 1 2 u inary, Euclidean) frequency of finite temperature T field 

d xdT(-(\7(j>) +-c {d r <t>) +-r<t> +-<f> ) (B2) thcory for fcrmion given by Um = (2m + 1)ttT and t x ^' 



1\ f 

5*0 = A / d 2 X / dT(f) *r.,aT**n,o 
J J n=l,o=l 

Here, the 2x1 Nambu spinor is given by 



*'n,a( k .<*>m) = (c+«,a( k , w m)> C_ n , a (-k, -W m )) 



are the usual Pauli spin matrices and r* is 2 x 2 identity 
matrix. 

(B3) Let us now however consider the zero temperature T = 

version of the same field theory in which the we convert 
the discrete Euclidean imaginary Matsubara frequency 
co m into continuous Minkowskian real frequency uj = fc 
by Wick rotation ioj m — > uj and the corresponding sum 
into integral. We now compute the one-loop Feynman 
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diagram for fermion self energy given by the following FIG. Bl One-loop Feynman diagrams for fermion self energy 
Feynman diagrams; and Yukawa vertex 




with the corresponding expressions, 



^lloop, % _ 1 f d3k i(Po + k )r z + v F {pi + k^T 1 - i](p 2 + k 2 ) 2 ir v 



* [P) ~NjJ (2tt)3 [(po + fc )2 + ^ + k,y + rf(p 2 + fc 2 )4](r + T 2 (fc)) 1 J 



where r 2 (fc) is the two-point function given in Eq. (8). 



r (k] = 4 / rf3 P + ") + ^fefe + fcx) + C 2 p>^ + fc,) 2 , A , , 

T ' ' (2tt)3 ((ft + w)a + t&(p B + fc x ) 2 + eiPv + k v mn 2 + vlpl + rfpf) + ;/ ' 



This is the fermion self energy correction from its cou- late the i(po + ko)r z + vf(.Pi + fci)r* part and try to 

pling to bosonic Ising nematic scalar field. Note that consider it as a 1 + 1 Dirac action. The Yukawa vertex 

we have not only the non-Lorentz p\ term but also the at zero external momenta-frequency is given by 
Lorentz symmetry breaking vp ^ 1 term even if we iso- 

m = W f j ^ rr [ TtG *^'P>'") TtG *^^' n ) Tt r + r a 1 (p,n) ) (B6) 



i 

These two quantities were used to analyze the structure Appendix C: Details on Quasiparticle Spectral 

of RG equations of the theory (equations 1, 2 and 3) in Function Calculation 

Section III. 

Here we give the expression for the elements of quasi- 
particle self energy needed for the calculation of its spec- 
tral function in in Section IV. From 



v . , , v - / d 3 P i& + ^)7o - VF(k x +Pxhi - y{k y + P y ) 2 l2 (ru 



where r\ = and A is the fermion-nematic coupling critical fluctuations on quasiparticle spectral function in 
constant in Eq.(3) which we recover here (rather than Section IV. Using representation (4), we have 
normalized to 1) in order to study the effect of nematic 



[ x ' y ' ' 1 (2^[{u J + ny + vi(k x +p x Y + r 1 2 {k v +p v Y}{T + T 2 {p)) ^ 
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z b {k x ,k y ,uj) = y 



d 3 p 



-v F (k x +Vx) 



(2tt) 3 [(w + ft) 2 + v 2 F (k x + Px y + r) 2 (k y +p y ) 4 ](r + T 2 (p)) 



(C3) 



E c (k x , k y ,Lo) 



d 3 p 



-V(ky+Py) 2 



(2tt) 3 [(W + ft) 2 + «2 (fe s + 7^ + p y )4]( r + r 2 (p)) 



(C4) 



These expressions are used to obtain the quasiparticlc 
spectral function shown in Fig. 7. 
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